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When heavy Quarkonia are formed in collisions between between nu¬ 
clei, their production cross section is modified relative to that inp + p col¬ 
lisions. The physical effects that cause this modification fall into two cat¬ 
egories. Hot matter effects are due to the large energy density generated 
in the nuclear collision, which disrupts the formation of the quarkonium 
state. Cold nuclear matter effects are due to the fact that the quarkonium 
state is created in a nuclear target. I will review experimental aspects of 
quarkonia production due to both hot and cold matter effects. 
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1 Introduction 


Interest in studying qnarkonia production in high energy heavy ion collisions was 
motivated decades ago by the prediction that J/^fj formation would be suppressed by 
color screening effects in the Quark Gluon Plasma (QGP) [I]. But experience has 
shown that relating J/^ suppression to the energy densities of the hot matter formed 
in heavy ion collisions is complicated by the presence of competing effects that also 
modify J/^fj production. The importance of this is illustrated in Figure where the 
nuclear modification factor, Raa-) is plotted versus collision centrality for collisions 
of heavy nuclei at ^/s^ = 17.3 to 200 GeV. There is no obvious dependence on the 
energy density produced in the collision. The measured nuclear modification factors 
are similar at mid rapidity at the two widely different collision energies, but at 200 
GeV show more suppression at forward rapidity, where the energy density is slightly 
smaller than at mid rapidity. Interpreting the data requires a detailed understanding 
of the mechanisms that compete with color screening. 
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Figure 1: The nuclear modification factor versus collision centrality for collisions of 
heavy nuclei at energies of ~ (In+In, Pb+Pb), 19.4 GeV (S+U), and 200 

GeV (Au+Au). 

In general, quarkonia production is modified by both cold nuclear matter (GNM) 
effects and hot matter effects [2] GNM effects modify the yield or kinematic distri¬ 
butions oi J/'ll) produced in a nuclear target in the absence of a QGP. They include 
modification of parton densities in a nucleus E 111 El [g [T], breakup of the charmo- 
nium or bottomonium or its precursor QQ state in the nuclear target due to collisions 
with nucleons |Hl Ellin], transverse momentum broadening as the QQ traverses the 
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cold nucleus, and initial state parton energy loss m- CNM effects are expected 
to be strongly dependent on collision system, collision energy, rapidity and collision 
centrality. 

In hot matter, there is not only suppression of quarkonia yields by color screening 
effects that reduce the binding energy of the quarkonium state, but also enhancement 
of quarkonia yields due to coalescence of QQ pairs that are initially unbound, but 
which become bound due to interactions with the medium (see for example [12] )• In 
particular, at LHC collision energies the yield of J/'ip from coalescence of a c and c 
from different hard processes becomes dominant. 

When discussing modihcation of quarkonia yields in nucleir collisions, the hi¬ 
erarchy of relevant time scales is very important. The nuclear crossing time for 
Lorentz contracted nuclei at RHIC(LHC) energies is ~ 0.3(0.001) fm/c. The forma¬ 
tion time of a J/ip meson in its rest frame is ~ 0.3 fm/c. The QGP thermalization 
time is ~ 0.3 — 0.6 fm/c. The QGP lifetime in central collisions at RHIC energy is 
~ 5 — 7 fm/c, longer at LHC energies. The J/'ip lifetime is ~ 2000 fm/c. Thus the 
creation of quarkonia and their modification in the hot medium occur on different 
time scales. They are often taken as factorable, although there are some theoretical 
models in which they can not be factorized m- 

2 proton-Nucleus collisions as a reference for CNM 
effects 

CNM effects are observed in either p-|-A or d-\-A collisions, which allow us to study 
the modihcation of quarkonia yields when they are produced in a nuclear target. It 
should be noted here that p-|-A collisions involving heavy targets appear to produce 
a small hot spot containing a large energy density, as evidenced by the observation 
of collective how effects m- However when studying hard probes (jets or heavy 
quarkonia), p-|-A collisions differ from A-l-A collisions. In p-|-A collisions the hot 
matter and the hard probes are both produced in a transverse area dehned by the 
radius of the projectile proton. In A-l-A collisions however a high energy density is 
produced throughout the collision volume, and the hard probe is typically immersed 
in that extended volume of evolving hot matter for several fm/c or more. We will 
return later to the question of whether the small hot spot produced in p-|-A collisions 
signihcantly modihes the yield of hard probes. 

There have been several studies of p-|-A or d+A data in which px integrated J/%p 
Raa values were corrected for shadowing effects using a shadowing parameterization, 
and then the data were fitted to extract an effective absorption cross section, aabs, 
where aabs accounts for all effects that modify J/ip production aside from shadowing. 
If it is assumed that CNM effects can be factorized from hot matter effects, then 
heavy ion collision data can be corrected for CNM effects using this parameterization 
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of the p+A data. 

In one such study (see section 5 of m) d+Au Rcp data from PHENIX at = 

200 GeV were corrected for shadowing effects using the EPS09 parameterization and 
values of aabs extracted at rapidities of —2.2 < y < —l.2y, y < 0.35 and 1.2 < p < 2.2. 
The parametrization was then used in a Glauber model to estimate the contribution to 
the Au+Au Raa from GNM effects at p < 0.35 and 1.2 < \y\ < 2.2, and the measured 
Au+Au Raa was divided by the GNM Raa- The result is shown in Eigurej^ (left). 
Comparison with Eig. shows that the correction for GNM effects eliminates the 
large difference in suppression for the Au+Au data at mid and forward rapidity. 
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Figure 2: Left: The measured Raa for J/tp production at 200 GeV at p = 0 and 
y = 1.7 divided by the estimated Raa due to CNM effects. Right: Comparison of the 
measured Raa divided by the estimated Raa due to CNM effects for J/^ production 
in Pb+Pb collisions at sqrtsn = 17.3 GeV at p = 0.5 and in Au+Au collisions at 
■sjspppp — 200 GeV at p = 0 

In a similar study of p+Pb data from the NA60 experiment at = 17.3 GeV, 

the Raa was corrected for shadowing effects using the EKS98 parameterization, and 
then the value of aahs extracted [H]. This was used to estimate the contribution to 
the Raa in Pb+Pb collisions so that the measured Pb+Pb Raa could be corrected 
for CNM effects. In Figure (right) the resulting CNM corrected results at 
= 17.3 GeV and y = 0.5 are compared with the CNM corrected mid rapidity data 
for Au+Au at = 200 GeV. The data are presented as a function of dN/dy, 

which is used as a proxy for energy density. It is found that when the A+A data 
are corrected for CNM effects the remaining nuclear modification seems to scale with 
energy density, and is about 50% for central collisions at RHIC energy. 
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3 Nucleus-Nucleus collisions and coalescence 


As mentioned earlier, in addition to suppression of the quarkonia due to color screen¬ 
ing, there is another effect that can occur due to QGP formation - coalescence. There 
are two coalescence scenarios. 

In the first, a Q and Q that were produced in the same hard process, and are nearly 
bound, can become bound through interactions with the medium. This first scenario 
occurs even if only one heavy quark pair is produced in a collision, and is present to 
some degree at all collision energies. It is sometimes referred to as ’’regeneration”. 

In the second scenario, a Q and Q that were produced in different hard processes 
can thermalize in the medium and combine statistically at hadronization. This second 
scenario becomes important only if a single collision produces many heavy quark 
pairs. At LHC energies we expect that ~ 100 cc pairs are created in a central Pb+Pb 
collision. 

Figure 1^ shows a comparison of measured Raa values for = 2.76 TeV ALICE 

Pb-|-Pb results and = 200 CeV PHENIX Au-|-Au results. The data, measured 

at forward rapidity in both cases, are plotted versus dNch/dr], which is used as a 
proxy for energy density. The suppression for the lower energy data is far stronger 
than for the higher energy data, suggesting that coalescence is now dominant. Indeed, 
the increase in J/tf) yield at the higher energy is concentrated at lower transverse 
momentum, as predicted by models that include coalescence at LHC energies [T5] . 


forward rapidity 



Figure 3: Comparison of Raa data at = 2.76 TeV by ALICE with 

Raa data by PHENIX at 200 CeV. 
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Of the A+A data available so far, the strongest suppression occurs at RHIC energy, 
presumably because at the higher LHC energies the gain in J/ip yield overcomes the 
loss due to stronger color screening effects. It is interesting to look at the dependence 
of the suppression for Au+Au on collision energy. Figure]^ shows Raa vs centrality for 
Au+Au collisions at ^/spppp = 39, 62 and 200 GeV [16]. The data are compared with 
theoretical calculations that include coalescence [I2|. The measured Raa is found to 
be similar at all three energies, and in the calculations that similarity is expected 
as a result of the stronger suppression due to color screening being balanced by an 
increased coalescence component as the energy increases. 

In addition to the Au+Au data, there are preliminary U+U J/ip data at = 

193 GeV from PHENIX at forward rapidity. They show slightly weaker suppression 
for central collisions than is seen in Au+Au data. This is consistent with a picture in 
which the increased suppression from the higher energy density in the U+U collision 
is more than balanced by the increased charm production from the larger number of 
nucleon-nucleon collisions in U+U n- 
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Figure 4: The Raa for Au+Au collisions at = 39, 62 and 200 GeV from 

PHENIX, compared with calculations that include both suppression due to color 
screening and coalescence. 


4 CNM effects in proton-Nucleus collisions 

There was a comprehensive study of fixed target J /ip p+A data from five fixed target 
experiments [ H | in which the measured Raa values for centrality integrated data were 
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corrected for shadowing effects using the EKS98 parameterization, and then fitted 
with an effective absorption cross section. Additionally there is a value of aabs for the 
J/^ from EKS98 corrected p+Pb data at = 17.3 GeV from [S]. A similar study 

oi J/'tp Raa values from y/s~^ = 200 GeV d+Au collisions at nine rapidities, this time 
corrected for shadowing using the EPS09 parameterization, was made in [10]. The 
aabs values extracted from all seven data sets were plotted versus the nuclear crossing 
time - the proper time in the rest frame of the evolving charmonium state during 
which it is exposed to collisions with nuclei. The results from PI are presented in 
Figure]^ 



Figure 5: Effective absorption cross sections extracted from shadowing corrected p+A 
or d+A data at a range of collision energies from = 17.3 to 200 GeV, plotted 

versus proper time spent in the target nucleus by the evolving charmonium state. 
The curve is discussed in the text. 

For nuclear crossing times of r > 0.05 fm/c the data at different energies seem 
to scale well with r. The data at r > 0.05/m/c were fitted [IHI with a model of an 
expanding color neutral meson passing through the target [iH] , showing that the data 
in this r range are consistent with breakup of charmonia. For shorter nuclear crossing 
times no such scaling is observed. This is expected, since those nuclear crossing times 
are shorter than the cc formation time, suggesting that for short nuclear crossing times 
the effective absorption cross section reflects physical processes other than breakup 
of the charmonia, such as energy loss in cold nuclear matter. 

The observation of elliptic flow in p+A collisions at LHG energies and d+A colli¬ 
sions at RHIC energies (see for example PI) suggests that a small fireball is produced 
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in those collisions. This raises the question of whether quarkonia production in p+A 
collisions is modified by final state effects. 

For the the p+A data from ALICE at the LHC seem to be consistent with 
models of CNM effects, including coherent energy loss at forward rapidity [IH], al¬ 
though there are significant experimental and theoretical uncertainties. At RHIC 
energy, with similar caveats, there is no strong evidence that more than CNM effects 
are needed to describe the PHENIX d+Au data [11]. It is worth noting also that 
the scaling observed with r in Eig. |^for r > 0.5 fm/c would be broken if there were 
strong final state effects in d+Au collisions at 200 GeV, but it seems to hold within 
the uncertainties on the extracted aabs values. In all cases, however, because feed 
down from the 'tp' accounts for only about 10% of the inclusive J/ip yield, it is still 
possible that the very weakly bound xp' is significantly modified by final state effects. 

Unexpectedly strong suppression of xp' yields in central d+Au collisions was ob¬ 
served at ^/spppp = 200 GeV at mid rapidity by PHENIX (20], as shown in Figure A 
very similar strong suppression of the xp' yield is observed in central p-l-Pb collisions 
at y/spppp = 2.76 TeV at both forward and backward rapidity in ALICE preliminary 
data EH. The strong differential suppression of the xp' relative to the J/xp can not be 
due to a larger nucleon breakup cross section for the larger, more weakly bound xp' 
because the nuclear crossing time is too short at all of the collision energies and ra¬ 
pidities of the data. Nor is it expected from any other known CNM process. It seems 
likely that the strong xp' suppression is caused by final state effects. Destruction of the 
xp' by co-moving particles in the final state has been proposed as an explanation [22] . 


5 Upsilons 

The comparison of charmonium measurements at RHIC and LHC energies has re¬ 
vealed some very interesting physics. However it has not provided a direct comparison 
of color screening effects at RHIC and LHC temperatures because the dominant mech¬ 
anism for charmonium production at LHC energies, charm coalescence, is different 
from that at RHIC energies, where color screening dominates. 

The study of Upsilon production in A-l-A collisions offers several advantages over 
charmonium. Eirst, the T(1S), T(2S) and T(3S) states all have significant branches 
to dileptons, and can thus be observed in the same experiment. The three states 
span a broad range of binding energies and sizes, and have similar shadowing effects. 
Importantly, the T states will not have large coalescence contributions at RHIC or 
LHC because the bottom production rate in central collisions at LHC energies is 
similar to the charm pair production rate at RHIC. On the down side, the T pro¬ 
duction cross sections are small, and the difference in mass between the states is 
also small. Therefore mass resolved T measurements require large luminosity, large 
detector acceptance, and excellent momentum resolution. 
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Figure 6: The RdAu for 'tp' production in d+Au collisions at = 200 GeV mea¬ 

sured at mid rapidity by PHENIX. The RdAu for J/i’ production is shown for com¬ 
parison. 

CMS has measured the ratios of T(2S) and T(3S) to T(1S) in p-|-p and p-l-Pb 
collisions [23]. In p-|-Pb collisions they observe that the T(2S) and T(3S) are differ¬ 
entially suppressed relative to the T(1S) by about 17% and 30% respectively. The 
differential suppression of the excited states is found to be stronger for events with 
larger particle production, suggesting that the larger number of particles in the final 
state has a stronger effect on the more weakly bound states. 

The CMS experiment has measured the T modification in Pb-|-Pb collisions at 
= 2.76 GeV [21|. The Raa values are shown in Figurej^for the T(1S) and T(2S) 
states. The T(3S) state is so strongly suppressed that values were not extracted. The 
fact that the (2S) and (3S) states are much more strongly suppressed than the (IS) 
state is suggestive of stronger color screening suppression for the larger, more weakly 
bound excited states. The data are well described by models in which an T potential 
model is modified by color screening in an evolving quark gluon plasma [2Sl[2n|. The 
statistical precision of the LHC data is expected to improve greatly during future 
planned running of the LHC program. 

The measurement of Upsilons at RHIC energy is very challenging because the 
production cross sections are small. So far, the measurements have been made with 
low statistical precision. The three states are unresolved in PHENIX, while in STAR 
only the T(1S) state can be separated by line-shape fitting. Therefore the existing 
RHIC measurements [271 EH] do not provide strong constraints on theoretical models. 
This will change in the future with data from the new STAR MTD detector [29], 















which was operational in the 2014 RHIC rnn, and future data from the proposed 
sPHENIX experiment at RHIC [SO]- 
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Figure 7: The Rppt for T production at = 2.76 TeV measured at mid rapidity 

by CMS. The theory curves are from [26]. 

References 

[1] T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 

[2] N. Brambilla et aL, Eur. Phys. J. C 71, 1534 (2011) 

[3] M. Hirai, S. Kumano and T.-H. Nagai, Phys. Rev. C 76, 065207 (2007) 

[4] K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0904, 065 (2009) 

[5] K. Kovarik, I. Schienbein, E. I. Olness, J. Y. Yu, C. Keppel, J. G. Morfin, 
J. F. Owens and T. Stavreva, Phys. Rev. Lett. 106, 122301 (2011) 

[6] D. de Florian, R. Sassot, P. Zurita and M. Stratmann, Phys. Rev. D 85, 074028 

( 2012 ) 

[7] 1. Helenius, K. J. Eskola, H. Honkanen and C. A. Salgado, JHEP 1207, 073 

( 2012 ) 

[8] C. Lourenco, R. Vogt and H. K. Woehri, JHEP 0902, 014 (2009) 


9 













[9] R. Arnaldi et al. [NA60 Collaboration], Phys. Lett. B 706, 263 (2012) 

[10] D. C. McGlinchey, A. D. Frawley and R. Vogt, Phys. Rev. C 87, 054910 (2013) 

[11] F. Arleo and S. Peigne, Phys. Rev. Lett. 109, 122301 (2012) 

[12] X. Zhao and R. Rapp, Phys. Rev. C 82, 064905 (2010) 

[13] B. Z. Kopeliovich, 1. K. Potashnikova, H. J. Pirner and 1. Schmidt, Phys. Rev. 
C 83, 014912 (2011) 

[14] V. Khachatryan et al [CMS Collaboration], Phys. Rev. Lett. 115, 012301 (2015) 

[15] J. Adam et al [ALICE Collaboration], arXiv: 1506.08808 [nucl-ex]. 

[16] A. Adare et al [PHENIX Collaboration], Phys. Rev. C 86, 064901 (2012) 

[17] D. Kikola, G. Odyniec and R. Vogt, Phys. Rev. C 84, 054907 (2011) 

[18] E. Arleo, P. B. Gossiaux, T. Gonsset and J. Aichelin, Phys. Rev. G 61, 054906 

( 2000 ) 

[19] A. Adare et al [PHENIX Collaboration], Phys. Rev. C 87, 034904 (2013) 

[20] A. Adare et al [PHENIX Collaboration], Phys. Rev. Lett. Ill, 202301 (2013) 

[21] R. Arnaldi [ALICE Collaboration], Nucl. Phys. A 931, 628 (2014) 

[22] E. G. Eerreiro, Phys. Lett. B 749, 98 (2015) 

[23] S. Ghatrchyan et al [GMS Gollaboration], JHEP 1404, 103 (2014) 

[24] S. Ghatrchyan et al [GMS Gollaboration], Phys. Rev. Lett. 109, 222301 (2012) 

[25] A. Emerick, X. Zhao and R. Rapp, Enr. Phys. J. A 48, 72 (2012) 

[26] M. Strickland and D. Bazow, Nncl. Phys. A 879, 25 (2012) 

[27] A. Adare et al [PHENIX Collaboration], Phys. Rev. C 91, 024913 (2015) 

[28] R. Vertesi [for the STAR Collaboration], arXiv: 1509.05359 [hep-ex]. 

[29] L. Ruan et al, J. Phys. G 36, 095001 (2009) 

[30] A. Adare et al, arXiv: 1501.06197 [nncl-ex]. 

[31] J. Adam et al [ALICE Collaboration], arXiv: 1506.08804 [nucl-ex]. 

[32] B. B. Abelev et al [ALICE Collaboration], Phys. Lett. B 740, 105 (2015) 


10 


